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Heavy traffic congestion daily occurs at merging sections on highway. For releasing this congestion, 
possibility of alternative configuration of vehicles on multiple-lane road is discussed in this paper. 
This is the configuration where no vehicles move aside on the other lane, ft has a merit in making 
smooth merging at an intersection or a junction due to so-called the "zipper effect". We show, 
by developing a cellular automaton model for multiple lanes, that this configuration is achieved by 
simple local interactions between vehicles neighboring each other. The degree of the alternative 
configuration in terms of the spatial increase of parallel driving length is studied by using both 
numerical simulations and mean field theory. We successfully construct a theoretical method for 
calculating this degree of the alternative configuration by using cluster approximation. It is shown 
that the theoretical results coincide with those of the simulations very well. 



I. INTRODUCTION 

In recent years, traffic dynamics has attracted much 
interest of pfiysicists, and thus it has been studied more 
and more diligently [l| 0. Researchers have mainly de- 
veloped the analysis of the traffic flow on one-lane road 
by using continuous models [3] and cellular automaton 
(CA) models [4]. Recently, the analysis of the traffic flow 
on multiple-lane road at an intersection or a junction is 
considered as one of the important study of traffic flow 
for releasing traffic congestion. Kita modeled merging 
interactions with game theory ^ . Hidas investigated ve- 
hicle interactions in merging and weaving by using agent 
based simulations [bI]. Davis introduced the cooperation 
in merging by adding interactions between pairs of ve- 
hicles in opposite lanes He showed that velocity of 
vehicles in cooperative merging was higher than that in 
no cooperation. 

However, these previous works were qualitative and 
did not study the configuration of vehicles in detail on 
two lanes before an intersection or a junction which de- 
termines the efficiency of merging. Among various con- 
figurations, the alternative configuration is the best be- 
cause it realizes the smooth "zipper" merging, which is 
the merging of vehicles by turns on one lane and on the 
opposite lane. Thus, the transformation of the configu- 
ration toward the alternative configuration is significant 
for the improvement of the congested flow at a junction. 

The purpose of this paper is to propose a simple and 
natural method for achievement of the alternative con- 
figuration of vehicles on two lanes. Then we study the 
method by both computer simulations and a mean field 
theory. The communication of vehicles on two lanes is 
studied by the CA model which we call the multiple lanes 
stochastic optimal velocity (MLSOV) model. This is an 
extension of SOV model proposed in Q by introducing 
the interactions between vehicles on the opposite lane. 



This paper is organized as follows. In Sec. II we define 
MLSOV model, and results of simulations of two-lane 
flow before an intersection is presented in Sec. IH. Then 
in Sec. IV we calculate the degree of alternative config- 
uration by using cluster approximation and compare it 
with the results of the simulations. Sec. V is devoted to 
the concluding discussions. 



II. A MULTIPLE LANES CA MODEL 
WITHOUT LANE CHANGE 

We propose here the MLSOV model by introducing in- 
teractions between two lanes. We choose the SOV model 
as a basis of our extended model because it is one of the 
simplest model whose fundamental diagram, i.e. fiow 
versus density plot, shows the metastable state observed 
in real traffic data @- MLSOV model represents the in- 
teractions of vehicles on two lanes by seeing the ones on 
the opposite lane each other. Note that we do not take 
into account of the lane change behaviour in this paper 
because we focus only on the achievement of the alterna- 
tive configuration of vehicles before an intersection. In 
the model, the movement of all vehicles on both lane 1 
and lane 2 is ruled as follows, i-th vehicle on each lane 
moves one cell in front in one time step with probability 
w* at time t provided that the next cell is empty. Thus 
the movement is described as 

^t+i _ I a:^* + 1, with probability 
^ \x\, with probability 1 — w|, 

where x\ is the position of i-th vehicle at time t as shown 
in figure [T] is called intension and normalized as < 
vl <1. z+l-th vehicle is located at the cell which is Ax* ^ 
cells ahead of i-th vehicle on the same lane as shown in 
figure [TJ The length of one cell is set as 7.5 m. 
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In the MLSOV model, the time evolution of vj is de- 
termined not only by i-th and i + 1-th vehicle but also 
j + 1-th vehicle which is the closest vehicle to i-th vehi- 
cle on the neighboring lane, with the distance i cells 
ahead as shown in figure [T] v\ is given as 

-a{nAa:*H,Ax*,)-«*}, (2) 

where V{/^x\ ^, ^^x\ J is the two-lanes Optimal Velocity 
(OV) function v\ tends to approach V{Ax\ Accj J 
with the increase of the sensitivity parameter a (0 < a < 
1). We define for convenience Ax^ ^ = oo in the case that 
i -f 1-th car does not exist and Ax\ ^ = oo in the case that 
j -I- 1-th car does not exist. 

Now, we set the concrete form of V as 
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FIG. 1: The vehicles that affect the behaviour of i-th vehicle 
in MLSOV model, i -I- 1-th vehicle is at the cell which is Aa;i j 
cells ahead of i-th vehicle on the same lane, j -I- 1-th vehicle is 
the closest vehicle to the i-th vehicle with the distance Aa;|_i 
cells ahead on the neighboring lane. 




(3) 



Ax*^^ > 1 and b.x\ .^ > 2, 



which is simple enough to achieve the alternative config- 
uration as shown below. Here we consider special two 
cases of the time evolution of u* with the initial condi- 
tion of = p. In the case a — 0, MLSOV model corre- 
sponds to the single-lane Asymmetric Simple Exclusion 
Process (ASEP) y\, since each vehicle moves irrespective 
of neighboring ones and is given as 




A-tL = 
A-T*, > 1. 



(4) 



In the case a = 1, MLSOV model corresponds to a two- 
lane version of the Zero Range Process (ZRP) [l3|, and 
is given as 




Ax[ ^ > 1 and Ax^ i = 
> 1 and Axl .^ = 1 
Ax*i , > 1 and Ax^ ^ > 2. 



(5) 



Note that in the case a = 0, the model is exactly solvable 
and we can theoretically calculate the stationary state 
and hence the fundamental diagram. In the case a = 1, 
if we neglect the effect of the other lane, then we can also 
solve the model exactly. 



III. SIMULATIONS 

Our idea for making the traffic flow smooth is to draw 
the compartment line between the lanes as shown in fig- 
ure m The line prohibits vehicles from changing lanes, 
and this will be used for communicating vehicles between 
lanes and expected to induce the zipper effect at the end 
of the line. Then, the disordered lane change at the merg- 
ing area will be smooth and becomes safe. 



The part of two-lane road indicated in figure [2] (b) is 
partitioned into identical cells as shown in figure [5] (c). 
The cell size is fixed and one cell can accommodate at 
most one vehicle. The boundaries are open, and each ve- 
hicle is updated in parallel, and it enters in the leftmost 
cell on lane 1 or lane 2, and moves straight ahead with- 
out changing lanes, and goes out of the rightmost cell. 
Although the vehicles cannot change lanes, vehicles on 
one lane can interact with the ones on the opposite lane 
by looking at the behaviour of these vehicles. 

All the parameters are defined in figure [2] (c). The 
length of this section is d, and the leftmost cells are at 
X — and the rightmost cells are at x = d — 1. 

In the simulations we choose the worst injection con- 
dition for merging, i.e., a vehicle enters in a; = simul- 
taneously on lane 1 and lane 2 with the probability a 
as long as both leftmost cells are empty. This injection is 
considered to be the most severe configuration for attain- 
ing the alternative configuration at the end of the road 
where merging occurs. We will demonstrate in the follow- 
ing that the alternative configuration is surely achieved 
in all cases of the injection. The probability of going out 
of the rightmost cell of lane i is set as /3i (i = 1, 2). 

Next, we define two important quantities in the simu- 
lations. Firstly, we quantify the degree of the alternative 
configurations of vehicles named Geminity {Ge). Ge is a 
function of x, and Ge{k) denotes the degree of the alter- 
native pattern of vehicles at x = k. In the simulations, 
Ge{k) is calculated by counting the state of the four cells 
at X = {k, k + 1}. There are 10 kinds of state labelled as 
S{ti) [n — 1,2,. ..10) shown as figure [H The symmetry 
between lane 1 and lane 2 is taken into account for elimi- 
nating the similar state. We defined c{n)k {n = 1, 2, ...10) 
as the total number of each S{n) at x = {k,k + 1} ap- 
peared through M times of simulations under the same 
conditions. The period for measurement of each sim- 
ulation is set as Ti < i < T2 — 1. When there is at 
least one vehicle at x — k, the state of the four cells at 
x = {k,k + l} can be n = 3, 5, 6, 7, 8, 9, 10 in S{n). Only 
5(3) among them represents the perfect alternative state 
which is defined by figure SI Thus, Ge(fc) is given by 
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FIG. 2: (a) An example of traffic flow on a weaving sec- 
tion before using our plan of drawing a compartment line at 
the merging area. Disordered lane changes cause traffic con- 
gestion at the merging area, (b) An example of traffic flow 
after introducing our plan. The line prohibits vehicles from 
changing lanes, and is expected to make smooth lane changes 
by achieving alternative configurations. (c)A CA model of a 
two-lane road with the compartment line. Each vehicle enters 
in the cell at a: = on both lane 1 and lane 2 simultaneously 
with the probability a, and goes out of the cell at a; = d — 1 
with the probability Pi {i = 1,2). 



c{n)k as 

Ge{k) = c(3)fe/(c(3)fe + c(5)fe + c(6)fe 

+ c{7)k + c{S)k + c{9)k + c(lO)fe) . (6) 

Ge ranges from to 1. The large value of Ge{k) denotes 
that the alternative configuration of vehicles at a; = A: is 
highly achieved. 

Secondly, we define the mean intension, i.e., mean ve- 
locity on the two cells ata; = fc (0<A;<d— 1) denoted 
by v{k), which is given by 



m = 



Y^A/-l Y^Wj-l Y^Ta-l (r 
2^j=a 2^i=0 2^t=Ti '^i"x\,k 
Y^Ta-l r ■ 
2^j"=0 2^i=0 2^t=Ti "x*.k 



(7) 



where S^t is defined bs 1 ii x\ — k, and if a;* ^ k. 
This is calculated through M times of simulations under 
the same conditions, and Nj (0 < j < M — 1) is the 
total number of vehicles entering in the leftmost cells on 
both lane 1 and lane 2 in j-th simulation. The period 
for measurement of each simulation is the same as that 
of Ge. 
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FIG. 3: 10 kinds of the state labelled by S{n) in the four 
cells at X = {fc,fc-|-l}. The symmetry between lane 1 and 
lane 2 is taken into account to reduce the number of states. 
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FIG. 4: The perfect alternative configuration on the two-lane 
road. When at least one vehicle exists ai x = k, only 5(3) 
represents this perfect alternative configuration. 



We obtained Ge{x) versus x and v versus x by 
using numerical simulations. The conditions of the 
simulations are follows. Three sets of parameters in 
the OV function are chosen as p = 1 and (g, r) = 
{(0.99, 0.99), (0.8, 0.8), (0.5, 0.5)}. Moreover, five kinds 
of the value of the sensitivity parameter are chosen as 
a = {0,0.001,0.01,0.1,1}. The length of the road is set 
as d = 100. The probability of simultaneous injection on 
both lane 1 and lane 2 is given as a = 0.05. Vehicles en- 
ter in the leftmost cells with the initial intension p which 
appears in ((21). Pi {i = 1,2) is given as each vehicle's 
intension v on the cell at x = d — 1. We see that traffic 
fiow is always in a free flow under this a and (3i. The 
number of iterations of the simulations is M = 10, and 
the period for measurement is between Ti = 100000 and 
Ta = 200000. 

The results of the simulations are shown in figure jS] 
From figure [Hj we see that Ge{x) increases monotoni- 
cally from to 1 as x became larger in the cases of 
a — {0.001,0.01,0.1, 1}, and stays constant at regard- 
less of X in the cases of a = 0. The sharpness of the in- 
crease of Ge{x) becomes larger as a becomes larger, and 
as q{= r) becomes smaller from figure|n](a)-(c). v{x) has 
one minimum value in the case of a = {0.001, 0.01, 0.1, 1}, 
and stays constant at v{x) = 1 in the case of a = from 
figure |6| (d)-(f). The position x which gives the mini- 
mum value of V becomes smaller as a becomes larger, 
and as q{= r) becomes smaller. The minimum value of 
v{x) becomes smaller as a becomes larger, and as q{— r) 
becomes smaller. 

Figure |S| (a)-(c) clearly show the achievement of the 
alternative configuration toward the spatial axis. These 
figures suggest us how we realize the smooth merging 
at an intersection or a junction. The sufficient length 
of compartment line for communicating is significant for 
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the achievement of the alternative configurations. For 
instance, if we set the target Ge as Getar = 0.9 in the 
case of a = 0.1, p = 1, and q — 0.5, which are typical 
values of real traffic [11], then we achieve this value of 
Ge by drawing the line of 22 cell length which is equal to 
165 m. Note that this alternative configuration is formed 
only by the local interaction within the range of Ax2 < 1 
of each vehicle. The relationships between a and the 
sharpness of the increase of Ge{x) suggest that the length 
of line necessary for realizing the complete alternative 
configuration becomes shorter if vehicles respond to the 
ones on the opposite lane more quickly. The relationships 
between q and the sharpness of the increase of Ge{x) 
suggest that this length of line becomes shorter if vehicles 
decelerate more strongly in the case of Ax2 < 1. 

The spatial change of v{x) is caused by the interactions 
between the vehicles neighboring each other in two steps. 
Firstly, i-th vehicle enters into a; = with Ax2,i = and 
slows down due to the OV function V < 1 to apart from 
the ones on the neighboring lane. It keeps slowing down 
until Ax2^i > 2. Then, Ax2,i = 2 is attained, and it starts 
to accelerate with V — 1. These interactions generate 
the minimum value of v{x) at the point where Ax2 of the 
most vehicle becomes 2. Vehicles slow down more quickly 
as the response parameter a becomes larger, and q{= r) 
becomes smaller. Thus, x which give the minimum value 
of v(x) becomes smaller as a becomes larger, and as q 
becomes smaller. 



IV. ANALYSIS OF ALTERNATIVE 
CONFIGURATIONS WITH CLUSTER 
APPROXIMATION 

In this section, we theoretically study the achieve- 
ment of alternative configurations by using cluster ap- 
proximation. In the approximation, the two-lane road 
as shown in figure [2] is divided into the four cells at 
X — {k, fc + 1} (fc = 0, 1, . . . , c? — 2) which is denoted by 
Cfc. Then, the stationary state on Gk is calculated in 
the order of fc = 0, 1, . . . , d — 2. The state vector of Gk 
at time t is defined as H*^ = {n(l)*^, n(2)*., . . . ,n(10)^.}, 
where Ii{n)\, {1 <n < 10) is the probability of Gk having 
the state Sn at time t as shown in figure [31 The state 
transition of Gk from TL\ to H^^^ is defined by using the 
state transition matrix whose size is 10 x 10. This 
state transition is given as 



p^n 



k- 



The stationary state 11^ is given as the solution of 



(8) 



(9) 



= 1. 



with the normalized condition X]ra=i n(^) 
We construct P^ to solve ([9|) in the order of fc = 
0, 1, 2, . . . , d — 2. Pfc is given in the following by the 
boundary conditions (B.C.), and by the dynamics of ve- 
hicles. 



The B.C. of each P^ (fc = 0, 1, 2, . . . , d - 2) are given 
as shown in figure [51 The left B.C. of Pq is given strictly 
as vehicles entering in the leftmost cells simultaneously 
with probability a provided that both leftmost cells are 
empty. The right B.C. of Pq needs approximations be- 
cause the stationary state of Cfc (1 < fc < d — 2) is not 
yet calculated. The right B.C. of Pq is approximated as 
a pair of vehicles existing on both cells at x = 2 with 
probability a. d is approximated as d = a/(l -|-a) which 
is the expected density of the loss system. The left B.C. 
of P/c (1 < fc < d — 3) is the stationary state of Gk-i 
which is given strictly by T^'^_i- The right B.C. of Pfc 
(l<fc<d — 3)is the stationary state of Gk+i which is 
not yet calculated. This stationary state is approximated 
by n^-^, i.e., the stationary state at x = {fc -I- 1, fc -I- 2} 
is approximated as that at a; = {fc — 1, fc}. This approxi- 
mation represents well the spatial change of the configu- 
rations if the flow is in free state. The left B.C. of Pd-2 
is the stationary state of Gds which is given strictly by 
n^g. The right B.C. of Pd-2 is given strictly as vehi- 
cles go out of the rightmost cell on lane z (i = 1, 2) with 
probability Pi. 




FIG. 5: The B.C. of each Pfc (fc = 0, 1, 2, . . . , d - 2). The 
right B.C. of Pq is approximated as a pair of vehicles existing 
on both cells at a; = 2 with probability a. a is approximated 
as ce = a/{l + a) which is the expected density of the loss 
system. The right B.C. of Pfc (1 < fc < d — 3) is the stationary 
state of Ck+i which is approximated by II^i. 

The dynamics of vehicles in this cluster approximation 
is given as follows. We update the intension of vehi- 
cles spatially to represent the spatial change of the mean 
intension v observed in the simulation results (figure [6l 
(d)-(f)). To realize the spatial change of the intension, 
we define Ui^k (0 < fc < d — 2) as the intension of i-th ve- 
hicle on Gk- To give the specific form of Ui.fc, we use Vk 
defined as the intension common to the vehicles on Gk, 
and Vk (0 < fc < d — 2) defined as the stationary mean 
OV function of Gk- Vk is calculated by using 11^ as 

Vk =ipm)T + 9n(6)r + 2rn(7)r + rn(9)r )/ 
(n(3)r + n(5)r + n(6)r + 2n(7)r+ 
n(8)r + 2n(9)r + 2n(io)r), (lo) 



and Vk is updated spatially as 



Vk+i = (1 - a)vk + aVk- 



(11) 
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Ui^k is given by using Vk, and by the configuration on Ck 
as 

Ui^k = (1 - a)vk + aV{Axii, Ax2i) 



(1 - a)vk, 
(1 - a)vk + ar, 
(1 - a)vk + aq, 
{1 - a)vk +ap, 



Axh = 

Axii > 1 and A.X2i = 
Axii > 1 and A.X2i = 1 
Axii > 1 and Aa;2i > 2. 



(12) 



Now the concrete form of Pk is given as follows. We 
define Pk (1 < i,j < 10) as the element of Pfc at i-th 
row, and j-th column. Pk is the probability of the 
transition of the state from S{j) to S{i). For instance, 
Pk (1, 1) = 1 — a, and Pk (7, 1) — a. The other elements 
of Pq, and Pfc (1 < fc < d — 3), and Pd-2 obeying the B.C. 
of the cells at a; = d — 1, such that (3j = Ui ,i_i (j = 1, 2) 
for i-th vehicle, are straightforward and shown in '11]. 

The degree of the alternative configurations at a; = fc 
(0 < fc < d — 2) in the stationary state is defined as Ge^ 
which is given by using 11^ as 



Get = n(3)r/(n(3)r + n(5)r + n(6)r 
+ n(7)r + n(8)r + n(9)r + n(io)r 



(13) 



We measure Ge^ versus x, and versus x, and com- 
pare them with the simulation results. The conditions 
of the cluster approximation are given similarly to the 
simulations, d, p, q, r, a, and a are the same to those of 
the simulations. (3j {j = 1,2) is given as the Ui^d-i for 
each i-th vehicle. 

The results of the four cluster approximation are shown 
in figure [H] together with the simulation results to com- 
pare the theoretical results with the simulations. 

Both Ge!^ and Vx coincide with Ge{x), and v{x) re- 
spectively. The difference \Ge^ — Ge{x)\ is smaller in 
the case of a = and a — 1 than that of a = 0.1. 

The coincidence between Ge^ and Ge{x) with various 
a and q suggests that our four cluster approximation is 
a good theoretical approximation for obtaining the al- 
ternative formation. We can obtain the length of the 
communication line necessary for the target of Ge by 
the theoretical calculation as well as by the simulations. 
Various Ge^ are generated by the spatial update of the 
intension as shown in p0|) - p2| . \Ge^ — Ge{x)\ in the 
case of a = and a — 1 are smaller than that of a = 0.1 
because the calculation of Ui^k does not contain the ap- 
proximation of the spatial update of Vx in the case of 
a = and a = 1: Ui_k in the case of a = is given as 



and Ui.fc in the case of a = 1 is given as 



(14) 



0, Axi, = 

r, Axii > 1 and Ax2i = 

q, Axii > 1 and Ax2i — 1 

p, Axii > 1 and Ax2i > 2. 



(15) 



The coincidence of Vx and v{x) suggests that this spa- 
tial update of the intension is a good approximation for 
expressing the spatial change of the mean intension. 



V. CONCLUDING DISCUSSIONS 

In this paper we have proposed a simple method for 
the achievement of alternative configurations of vehicles 
by using MLSOV model. We obtain for the first time 
the achievement of alternative configurations of vehicles 
by emergent behaviour of drivers. This alternative con- 
figurations is achieved along a compartment line which 
is drawn between two lanes, and prohibits vehicles from 
changing lanes, and permits them to interact with the 
ones on the opposite lane. This achievement is significant 
for traffic engineering in realizing the smooth "zipper 
merging" , and is supported by the four cluster approxi- 
mation. Moreover, drawing compartment lines for real- 
izing this alternative configuration costs less than other 
methods for casing traffic congestion at an intersection, 
e.g., construction of a cubic interchange. Further im- 
provement of the flow at an intersection or a junction by 
the alternative conflguration will be published elsewhere 

El. 
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FIG. 6; Ge{x) versus x and v{x) versus x obtained by the simulations, together with Ge^ versus x and Vx versus x calculated 
by the four cluster approximation. Ge{x) and Ge'^ increase monotonically as x increases, and the sharpness of the increase of 
Ge{x) and Ge^ becomes larger as a becomes larger, and as q{=r) becomes smaller. v{x) and Vx have one minimum value. 
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